MENU 2007 

11th International Conference 

on Meson-Nucleon Physics and 

the Structure of the Nucleon 

September10-14, 2007 

IKP, Forschungzentrum Julich, Germany 

SCALAR RADIUS OF THE PION AND 
TWO PHOTONS INTO TWO PIONS. 
STRONG S-WAVE FINAL STATE 
INTERACTIONS 



t^^ , J. A. Oiler, L. Roca 

O I Departamento de Fisica 

pq ' Universidad de Murcia 

;> ■ Spain, E-30071 

O' 






C. Schat 

Departamento de Fisica, FCEyN, 
Universidad de Buenos Aires 
^ ■ Ciudad Universitaria, Pab.l, 

(1428) Buenos Aires, Argentina. 



Ph' 






(N 



O 



Abstract 



1^ . The quadratic pion scalar radius, (r^)J, plays an important role for present 

Psj ! precise determinations of vrvr scattering. The solution of the Muskhelishvili- 

Omnes equations for the non-strange null isospin (/) pion scalar form factor 

determines that (r^)^ = 0.61 ± 0.04 fm^. However, by using an Omnes 

representation of this form factor, Yndurain recently obtains (r^)^ = 0.75 ± 

t^ . 0.07 fm^. A large discrepancy between both values, given the precision, then 

results. We show that Yndurain's method is indeed compatible with the 

determinations from the Muskhelishvili-Omnes equations once a zero in the 

/\ ' scalar form factor for some S-wave 1 = T— matrices is considered. Once 

c^ ■ this is accounted for, the resulting value is (r^)J = 0.63 ± 0.05 fm^. 

On the other hand, we perform a theoretical study of the reaction 77 — > 
7r° based on dispersion relations. The large source of uncertainty for ^/s > 



TT 



0.5 GeV, due to variations in the phase used in the Omnes function above the 
KK threshold, is removed by taking one more subtraction in the dispersion 
relation. This allows us to make sharper predictions for the cross section so 
that one could use this reaction to distinguish between different low energy 
TTTT parameterizations, once independent experiments are available. We also 
study the role played by the a or /o(600) meson in this reaction and determine 
its width to two photons. 
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1 Introduction 

Here we summarize the two papers [1 , 2] that mainly handle with the strong 
influence of the 1 = S-wave meson-meson final state interactions. We 
concentrate here on the non-strange / = scalar form factor of the pion [1] 
and 77 -^ 7r*^7r*^ [2]. Both processes can be formulated in a way that has 
in common the same basic function in order to take care of the strong final 
state interactions in the / = S-wave. This function has been recently the 
origin of large uncertainties in its implementation in the literature, both for 
the scalar form factor of the pion [3-5] and for 77 — > vr^vr^ [6]. 

The scalar form factor of the pion, TT^{t), corresponds to the matrix ele- 
ment 

TM ^ /.Se-<.'-..>„')| K.(.).(.) . ™.3(.).(.)) W,)) , « ^ (,'-,)■■ 

(1) 
Performing a Taylor expansion around t = 0, 

r.(t) = r.(o)|i + it(r^): + o(t^)| , (2) 



where {r )^ is the quadratic scalar radius of the pion. The quantity {r )^ 
contributes around 10% to the values of the S-wave tttt scattering lengths 
a[] = 0.220 ± 0.005 M'^ and ag = -0.0444 ± 0.0010 M-\ as determined in 
Ref. [7] by solving the Roy equations with constraints from two loop Chiral 
Perturbation Theory (CHPT). If one takes into account that one has a pre- 
cision of 2.2% in the scattering lengths, a 10% of contribution from (r^)^ is a 
large one. Related to that, (r^)^ is also important in SU{2) x SU{2) CHPT 
since it gives the low energy constant £4 that controls the departure of F^^ 
from its value in the chiral limit [8, 9] at next-to-leading order. 

Based on one loop xPT, Gasser and Leutwyler [8] obtained (r^)J = 
0.55 ± 0.15 fm^. This calculation was improved later on by the same authors 
together with Donoghue [10], who solved the corresponding Muskhelishvili- 
Omnes equations with the coupled channels of tttt and KK. The update of 
this calculation, performed in Ref. [7], gives (r^)J = 0.61 ± 0.04 fm^. Mous- 
sallam [11] employs the same approach and obtains values in agreement with 
the previous result. One should notice that solutions of the Muskhelishvili- 
Omnes equations for the scalar form factor rely on non-measured T— matrix 
elements or on assumptions about which are the channels that matter. Other 
independent approaches are then most welcome. In this respect we quote the 
works [12-14] , and Yndurain's ones [3-5] . These latter works have challenged 
the previous value for (r^)J, shifting it to the larger (r^)^ = 0.75 ± 0.07 fm^. 
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If this is translated to the scattering lengths above, employing an equation 
of Ref. [7], it implies a shift of +0.006 Af-^ for a[] and -0.001 M'^ in ag. 
Thus, one is referring to a shift of slightly more than one sigma. Refs. [3,4] 
emphasize that one should have a precise knowledge of the / = S-wave 
phase shits, So{s), for s > 4M|. GeV^, M^ is the kaon mass, to disentangle 
which of the values, either that of Ref. [7] or [3], is the right one. However, 
this point is based on an unstable behaviour of the solution of Ref. [3] with 
respect to the value of 6o{4:Mj^). Once this instability is cured, as shown 
below, the resulting (r^)J only depends weakly on 5o(s), s > 4M|^, and is 
compatible with the value of Ref. [7] . 

Regarding the reaction 77 -^ 7r°7r° one has to emphasize that due to the 
absence of the Born term (as the vr'^ is neutral), this reaction is specially 
sensitive to final state interactions. For energies below 0.6 GeV or so, only 
the S-waves matter, which have J = or 2. It is in this point where both the 
study of this reaction and the scalar form factor match. Recently, Ref. [6] 
updated the dispersive approach of Ref. [15] to calculate (7(77 — > 7r°7r°). Here 
one finds a large uncertainty in the results for y/s > 0.5 GeV that at around 
0.6 GeV is already almost 200%. Again, this is due to the lack of a precise 
knowledge of the phase of the 77 ^ tttt / = S-wave amplitude above Arnj^. 

We showed in Refs. [1,2] that one can improve largely this situation by 
employing an appropriate Omnes function in the / = S-wave. The key 
point is that this function should be continuous under changes in the phase 
functions used above 1 GeV, a point overlooked in the previous studies. 



2 The scalar form factor 

Ref. [3] makes use of an Omnes representation for the pion scalar form factor. 



r,(() = P(()exp 



%.'- ^»(»'' 



^ J4M2 s'(s' -t-ie) 



(3) 



Here, P{t) is a polynomial in t normalized such that P{0) = Tt^{0) and whose 
zeroes are those of r^(t). On the other hand, 0o(^) is the continuous phase 
of r.,^{t)/P{t). Then Refs. [3,4] make use of asymptotic QCD which predicts 
that the scalar form factors should go as — 1/t times a positive smooth factor 
for t -^ +CXO, so that the phase of the form factor should tend to +tt in the 
same limit. At this point, Refs. [3,4] make an assumption that is not always 
necessarily fulfilled. Namely, to identify 0o(^) with the phase of r7r(t), that we 
denote in the following as p{t). If this identification is done, as in Refs. [3,4], 
it follows that P{t) must be a constant, r^(0), because the behaviour for 
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t — > +00 that follows from Eq. ([3]) is 

r^(t) ^ (-i)-*(°°)/"rt-<^(°°)/"r^(o) , (4) 

with n the degree of P{t). As QCD implies in this assumption that 0(oo)/7r = 
1, then n = and hence P{t) = Tt^IO), just a constant. One must be aware 
that in Eq. ([3]) 4>oit) is the phase of TT,it)/P(t). Notice that the phase of 
r7r(t) is not continuous when crossing a zero located at ti G M, as there is a 
flip in the sign when passing through. However, the phase of rT^{t)/P{t) is 
continuous, since the zero is removed. This is the phase one should use in the 
Omnes representation, Eq. (IHl), because it results from a dispersion relation of 
logr7r(t)/P(t), and then 0(t) must be continuous (but not necessarily p(t)). 
As stated, Ref. [3] took 



r,(t) = r,(o)exp 



t r ,,_ Pis') 



.^ Jam^ s'(s' — t — ie) 
So that the scalar form factor is given by. 



(5) 



The phase p{s) is flxed in Refs. [3,4] by invoking Watson's flnal state theo- 
rem. For s < sk, sk = 4M|^, it implies that p(s) = So{s), where neglecting 
inelasticity due to multipion states, an experimental fact. For 1.42 > ^/s > 
1.1 GeV, Ref. [3] stressed the interesting fact that experimentally the in- 
elasticity turns out to be small and hence Watson's flnal state theorem can 
be applied approximately again. In the narrow region between 2Mk and 
1.1 GeV inelasticity cannot be neglected but Ref. [3] argues that, as it is so 
narrow, its contribution to Eq. ([6]) is small anyhow and, furthermore, that the 
elasticity parameter rj is not so small, so that one could still apply Watson's 
flnal state theorem with corrections. Thus, for s^ < s < 2 GeV^, Ref. [3] 
identifles again p{s) ~ So{s). Finally, for s > sq = 2 GeV^ Ref. [3] takes a 
linear extrapolation from 5o(so) to tt. One should here criticize that it is still 
a long way to run from values of 5o(so) < 27r up to tt at s ^ +00. With all 
these ingredients, and some error estimates, the value (r^)^ = 0.75 ±0.07 fm^ 
results [3,4]. 

As discussed above in the lines of Ref. [1], the steps performed in Ref. [3] 
are not always compatible. In Ref. [1] we took as granted the assumption that 
Watson's flnal state theorem can be approximately applied for 1.5 GeV > 
^/s > 2Mk- Our assumption is in agreement with any explicit calculation of 
the pion non-strange 1 = scalar form factor [7,10,11,13] and it is the proper 
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generalized version of the assumption of Refs. [3,4] of identifying p(s) ~ 5q{s). 
Now, Watson's final state theorem implies that 0(s) = ip{s) (modulo vr), 
with ip{s) the phase of the / = S-wave tttt amplitude, t^^T, = (rje^''^^ — l)/2z. 
It occurs, as stressed in Refs. [4, 16], that (p{s) can be either ~ So{s) or 
~ So{s) — vr depending on whether 6o{sk) > tt or < tt, respectively, for 
sk < s < 2 GeV^. The latter case corresponds to the calculation in Ref. [7], 
while the former is the preferred one in Ref. [4] and arguments are put forward 
for this preference in this reference. 



Let us evolve continuously from one situation {6o{sk) < tt) to the other 
{Sq{sk) > tt). In the first case ip{s) has an abrupt drop for s > Sk simply 
because then rj < 1 and while the real part of 1^^^^ rapidly changes sign, its 
imaginary part is positive (> 0). The rapid movement in the real part is due 
to the swift one in 6o{s) in the KK threshold due to the /o(980) resonance. As 
a result for s < sk, ^{s) = So{s) ~ vr and for s > sk then ip{s) < 7i/2. This 
rapid movement gives rise to a rapid drop in the Omnes function, Eq. ([5]), so 
that the modulus of the form factor has a deep minimum around sk- Here, 
one is using Watson's final state theorem with (f)o{s) = (p{s) and the form 
factor of Ref. [10] is reproduced. Notice as well that in this case the function 
(j){s) approaches vr from below for asymptotic s and then P(t) = Tq{0) in 
Eq. ([3]). Now, we consider the limit So{s) -^ n^ for s -^ s^. The superscript 
— (+) indicates that the limit is approached from below(above). In the limit, 
the change in sign in the real part of t^r^ occurs precisely at sk, so that for s = 
s~^, (p{s) =11 and for s = s~^ then ip{s) < 7i/2 (indeed it can be shown from 
unitarity that must be 0). As a result one has a drop by — vr in ip{s) which 
gives rise to a zero in the Omnes representation of the scalar form factor. 
Thus, the deep has evolved to a zero when 6o{sk) — ^ tt^. Because of this 
zero the proper Omnes representation now involves a P(t) = r^(0)(l —t/sx) 
and 0(s) is no longer ip{s) but ~ ip{s) + tt ~ So{s) for 2.25 GeV^ > s > sk- 
This follows simply because 0(s) is continuous. Thus, we go into a new realm 
where (j){s) ~ So{s) and the degree of P(t) is 1, so that T.„(t) has a zero at the 
point Si where (5o(si) = it and Si < sk- Note that only at Si the imaginary 
part of r^(t) is zero and this fixes the position of the zero [1]. We should 
emphasize here that if one uses Eq. ^ with 0(s) ~ 5q{s), as in Refs. [3,4], 
then in the limit 5q{s) —* 7r+ for s ^ s~^ the Omnes representation would give 
rise to |r7r(sft:)| = oo, while in the previously discussed limit of 5o('S) — ^ vr^ 
for s ^ s]^ one has |r^(sii:)| = 0. This discontinuity was corrected in Ref. [1] 
and it is the benchmark for a jump by one unit in the degree of P{t), a 
discrete function, in Eq. ([H]). 
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Hence for 5{){sk) > tt one has to use 



r,(t) = r,(o) ( 1 - ^ ) exp 



Ms') 



— I ds ' ' ' — — 

^ Jam^ s'(s' -t-ie) 



(7) 



with 0(s) ~ 5q{s) for s < 2.25 GeV^. The uncertainties in this approximation 
for s > sk are discussed in Ref. [1] and included in the final error in (r^)^. 
The estimation is based in diagonalizing the J = S-wave S-matrix for 
s < 2.25 GeV^, so that two elastic channels can be singled out [4]. We also 
remark that now 0(s) for 5q{sk) > tt must tend to 2n asymptotically so 
as to match with the asymptotic behaviour of r^(t) as — 1/t. In this way 
we have now a very soft matching with asymptotic QCD since for s around 
2.25 GeV^, 5q{s) ~ 27r. This was not the case in Ref. [3,4]. Notice that from 
our work it follows that the precise knowledge of the asymptotic behaviour 
of the phase of the form factor is not relevant as 0(s) can tend either to 2tt 
{8o{.sk) > tt) or to n (5o(sx) < ^), and the results are very similar. 
Our final value is 

(r2)^ = 0.63 ± 0.05 fm^. (8) 

The error takes into account different tttt / = S-wave parameterizations, 
namely those of Refs. [7] and [17], the error in the application of Watson's 
final state theorem above 1 GeV and up to 1.5 GeV, and the uncertainties in 
(f){s) given by asymptotic QCD for s > 2.25 GeV^. This value is compatible 
with that of Ref. [7], (r^)^ = 0.61 ± 0.04 fm^, and also with (r^)^ = 0.64 ± 
0.06 fm^ of Ref. [13] calculated from Unitary CHPT. 



3 The 



77 



A' 



reaction 



In this section we briefly review Ref. [2]. This reference extended the ap- 
proach of Refs. [6, 15] so as to be less sensitive to the phase of the / = 
S-wave 77 — > tttt amplitude above sk- For this phase one has a similar situ- 
ation to that of the scalar form factor of the pion, it can be either ~ Sq{s) or 
~ 5o(s) — TT for 1 < s < 2.25 GeV^ [2,6]. In the approach of Ref. [6] this orig- 
inates an uncertainty that raises dramatically with energy above 0.5 GeV, 
such for y/s ~ 0.6 GeV it is already 200%. 

Let us denote by Fj{s) the S-wave / = 77 — ;► tttt amplitude. The 
approach of Ref. [6,15] is based on isolating the left hand cut contribution of 
Fj which is denoted by Lj. These authors also employ the Omnes function 



^lis) 



exp 



+00 



ds - 



)i[s- 



^ JiMi S'{s' 



(9) 
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where (f)i{s') is the phase of Fj{s). For / = 2 by the application of Watson's 
final state theorem one has that 02 (s) = 62(3). For / = and s < sk, 
(t>o{s) = 6o{s). In the interval 1.5 > y/s > 1.1 GeV, 0o = ^o (modulo vr) 
because inelasticity is small again, as already remarked. Similarly as in the 
scalar form factor one can have because of the onset of inelasticity above 
2Mk and up to 1.1 GeV, that 0o is given either by ~ (5o or ~ (5o — tt. 

Ref. [6] then performed a twice subtracted dispersion relation of the func- 
tion {Ff{s)—Lj{s))/Qi{s). An important point to realize is that the previous 
function has no left hand cut and that Fj/Qj has no right hand cut. Making 
use of the Low's theorem, which implies that Li{s) is given by the Born term 
Bi[s) for s — >■ 0, one is only left with two subtraction constants to be fixed. 
One of these constants can be fixed by requiring that the 77 —>■ tt^tt^ S-wave 
amplitude, F/v(s), has an Adler zero around M^. The other one was fixed 
in Ref. [6] by requiring that the 77 -^ tt+tt^ S-wave amplitude, Fc{s), tends 
to the Born term Bc{s) for s -^ up to O(s^). One has to say that Ref. [6] 
did not include axial vector exchanges which indeed give rise to a term that 
vanishes for s — *> only linearly in s. This gives rise to a difference in the 
cross section of around a 30% at y/s ~ 0.5 GeV. 

In order to better handle the ambiguities in (j)o{s) above 1 GeV, Ref. [2] 
only uses ^ois) of Eq. ©0 for 0o(s) ~ 60(3) - vr for s > 1 GeV^. For the 
case (f>o{s) ~ So{s) above 1 GeV Ref . [2] employs 



ilo{s) = (1 ) exp 



. I'",,, Ms') 



IT 



4M2 



(10) 



and then a twice dispersion relation of [Fq — Lq)/Qq is performed. It is 
important to realize, as stressed in Ref. [2], that because of the first order 
polynomial in front of the exponential in Eq. ( ITOl) . one indeed has a three 
times subtracted dispersion relation for [Fq — Lq)/Qq. Recall that the latter 
is the original function used in Refs. [6, 15]. 

Because of this extra subtraction one can reduce dramatically the sensi- 
tivity to the 00 (s) above 1 GeV. The conditions used to fix the at most three 
subtraction constants that appear in our scheme are: i) Fn{s) — > for s — > 
with the slope fixed by one loop CHPT [18] (with an uncertainty of around 
15%), ii) Fcis) -^ Bc\s) + 0(s) with the rest fixed by one loop CHPT (with 
the same 15% of estimated uncertainty). The third condition is an upper 
bound to the value of the resulting cross section in the /o(980) region so that 
it is smaller than 200 nb. Notice that its experimental value is smaller than 
40 nb and, hence, we take here a very conservative uncertainty. 



""^We already know about the lack of continuity of r2o(s) when 5t^{sk) crosses tt when 
taking 4>o{s) given by ip{s) as in the case of the scalar form factor. 
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Figure 1: Final results for the 77 — > vr vr cross section. Experimental data are 
from the Crystal Ball Coll. [19], scaled by 1/0.8, as | cos6\ < 0.8 is measured and 
S-wave dominates. The lighter band corresponds to Ref. [7] while the darker one to 
Ref. [17]. The dot-dot-dashed line results after removing the axial vector exchange 
contributions, as in Ref. [6] with (po{s) ~ So{s) — vr for s > 1 GeV^. The band 
along each line represents the theoretical uncertainty. The dotted line is the one 
loop xPT result [18] and the dot-dashed one the two loop calculation [20]. 



We show in Fig. [T]our results together with the experimental points from 
Ref. [19]. The darker band corresponds to employ Ref. [17] for 6o{s) below 
1 GeV and the lighter one to use Ref. [7] . One sees that now with more precise 
data one should be able to distinguish between different low energy 6o{s) 
parameterizations as the theoretical uncertainty is much reduced. The widths 
of the bands correspond to the uncertainties related to the So{s) and 62(8) 
parameterizations used, those in fixing the three subtraction constants and in 
employing Watson's final state theorem for s > 1 GeV^, and it also includes 
the uncertainty in the asymptotic 0/(s) employed. In the figure we also show 
with the dotted line the one loop CHPT result [18] and with the dash-dotted 
line the two loop one [20]. There is a clear improvement when going from 
one to two loops in CHPT, though to have a perfect agreement with our 
results some higher order corrections are still needed. Finally, the dash- 
double-dotted line corresponds to the result of Ref. [6] with 4>o{s) ~ So{s) — vr 
for s > 1 GeV^. Let us recall that Ref. [6] does not include axial vector 
exchanges. Were they included, the results of this reference would fall inside 
the bands shown by our results. 
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By analytical continuation on the complex plane one can determine the 
coupling of the a to 77, gay-r, and calculate the width to 77 of this resonance 

= (2.1 ±0.2) X 10"3 , 

[6] corresponds to this 



ff<777 



[2]. We then obtain for the ratio of couplings 

with (yfo-^TT the (7 coupling to two pions. The result of 
ratio being 20% bigger at (2.53 ± 0.09) x 10^'^. Half of this difference is due 
to the omission of the exchanges of axial vector resonances in [6], and the 
other half comes from improvements delivered by our extra subtraction and 
our slightly different inputs. As a result, using the same value for Iganirl as 
in [6], our resulting value for r(o" -^ 77) would be around a 40% smaller 
than that in [6]. Taking into account different choices of Igaiml we end with 
r{a ^ '-f'-f) in the interval 1.8 — 3 KeV. 



4 Conclusions 

We have shown that both Yndurain's method [3] and the solution of the 
Muskhelishvili-Omnes equations [7, 10] provide compatible results for the 
quadratic scalar radius of the pion. The origin of the discrepancy between 
Refs. [3] and [7] was due to overlooking a zero in the scalar form factor in 
the former reference. We finally obtain [1] (r^)^ = 0.63 ± 0.05 fm^ and £4 = 
4.5 ±0.3. These numbers are in good agreement with (r^)^ = 0.61 ±0.04 fm^ 
and £4 = 4.4 ± 0.2 of Ref. [7]. 

We have also studied the 77 -^ tt^tt" reaction for energies ^/s < 0.7 GeV, 
where S-waves dominate. We have extended the original approach of Ref. [6, 
15] by performing a three times subtracted dispersion relation [2], instead of 
the twice subtracted originally employed. The sensitivity of the results with 
respect to the phase of the J = 77 -^ vrvr S-wave above 4M|- is then largely 
reduced. A key point is to properly handle the contribution of the /o(980) 
resonance, at least at the level of the order of magnitude. Importantly, 
one can then use this reaction to distinguish between different low energy 
TTTT parameterizations once new data on (7(77 — > tt^tt^) are available. The 
r(o" — i> 77) width is estimated in the range 1.8 — 3 KeV [2]. 
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